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We studied the finite-size effects on the magnetic behavior of the quasi-one-dimensional spin S 
= j Heisenberg antiferromagnets Sr 2 Cu 03 , Sr2Cu0.99M0.01O3 (M = Zn and Ni), and SrCuC> 2 . 
Magnetic susceptibility data were analyzed to estimate the concentration of chain breaks due to 
extrinsically doped defects and/or due to slight oxygen off-stoichiometry. We show that the sus¬ 
ceptibility of Sr2Cu0.99Ni0.01O3 can be described by considering Ni 2+ as a scalar defect (S e // = 0 ) 
indicating that the Ni spin is screened. In Sr2Cu0.99Z1i0.01O3 susceptibility analysis yields a defect 
concentration smaller than the nominal value which is in good qualitative agreement with crystal 
growth experiments. Influence of doping on the low-temperature long-range spin ordered state is 
studied. In the compound SrCu02, consisting of zigzag S = | chains, the influence of spin frustra¬ 
tion on the magnetic ordering and the defect concentration determined from the susceptibility data 
is discussed. 

PACS numbers: 75.10.Jm, 75.30.Cr, 75.10.Pq, 81.10.Fq 

Keywords: One-dimensional Heisenberg antiferromagnet, Spin chain, finite-size effects, pseudogap, chain- 
breaks, defects 


I. INTRODUCTION 


One-dimensional (ID) spin S = 1 Heisenberg antifer¬ 
romagnetic (HAF) model is the cornerstone of our under¬ 
standing of quantum magnets. While the exact eigenval¬ 
ues and eigenstates of this model can in principle be ob¬ 
tained using the Bethe ansatz 1 , calculations of thermody¬ 
namic quantities at nonzero temperatures has remained 
an active area of research. Probably, the first real at¬ 
tempt to calculate the susceptibility x(T) and the specific 
heat C P (T) of the S = \ HAF chain was by Bonner and 
Fischer (BF) 2 , who obtained numerical results for chains 
containing up to 11 spins subject to a periodic boundary 
condition and for temperatures T > HAJ/ks (where J is 
the intrachain nearest neighbor exchange and is the 
Boltzmann constant). These results revealed a broad sus¬ 
ceptibility maximum near T = 0.675 J/ks- At lower tem¬ 
peratures, the susceptibility decreases below the maxi¬ 
mum value due to the gradual appearance of large anti¬ 
ferromagnetic domains characterized by domain-wall-like 
S = \ excitations called spinons. Due to quantum fluc¬ 
tuations, the domain-wall dynamics persist down to zero 
temperature. Therefore, the ground state of the S = \ 
HAF chain is characterized as a spin-liquid in which the 
spin correlations decay as a power law 3 and the zero- 
temperature susceptibility x(0) has a finite value, which 
can be calculated exactly using the Bethe ansatz 4 . The 
%(T) behavior at low temperatures, where the BF re¬ 
sults are not accurate, was investigated by Eggert, Af¬ 
fleck, and Takahashi (EAT) using numerical Bethe ansatz 
(T > 0.003J/fcs) and field-theory methods 5 . The result¬ 
ing x(T) showed a surprisingly rich behavior with the 
existence of an inflection point near T = 0.087J/fcs, be¬ 
low which the slope of x(T) increases rapidly approaching 
infinity as T goes to zero. The field theory revealed a log¬ 


arithmic leading order correction to the zero-temperature 
susceptibility x(0). The accuracy of this result was fur¬ 
ther improved in the subsequent works by Lukyanov us¬ 
ing field theory 6 , and Kliimper and Johnston using the 
numerical Bethe ansatz 7,8 . 

Experimental confirmation of this logarithmic behav¬ 
ior is limited, because in most quasi-ID systems either 
the intrachain interaction is too small, making the tem¬ 
perature range of interest difficult to access, or the sys¬ 
tem undergoes a dimensional crossover at low tempera¬ 
tures due to interchain interactions (J’) (see, for example, 
Table 3 in ref. 9). Even in those rare cases where an ex¬ 
cellent one-dimensionality and exceptionally large J can 
be combined into a single material, the inevitable pres¬ 
ence of chain-breaks vitiate the ideal chain behavior at 
low temperatures. Due to chain breaks a chain is cut 
into finite-length segments. The staggered moment at 
the boundaries of these segments contributes a Curie-like 
susceptibility with an effective paramagnetic Curie con¬ 
stant that depends on temperature 10 . It has been argued 
that in a nominally pure crystal additional paramagnetic 
impurities play a minor role in giving rise to the Curie 
tail 11 . That the boundary effects are significant is also 
apparent from the fact that the measured susceptibility of 
powdered samples of quasi-ID chains can be significantly 
larger than that of the single crystals (for comparison see, 
for example, the susceptibility data published in Refs. 12 
and 13 on powder samples and that in Ref. 14 on single 
crystals of the same quasi-ID compounds). Therefore, 
in analyzing the experimental data of a nominally pure 
spin chain compound or of one doped with a dilute level 
of non-magnetic impurities within the chains, it is impor¬ 
tant to take the boundary contribution into account for 
an accurate description of the experimental susceptibil¬ 
ity. Conversely, knowledge of the boundary susceptibility 
allows for an accurate estimate of the defect concentra- 
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tion in a spin chain, which can otherwise be substan¬ 
tially underestimated 11 . Since any real systemmust have 
boundaries due to defects (or chain breaks), it is perti¬ 
nent to investigate and understand the effect of bound¬ 
aries on the physical properties and to test the validity of 
the theoretical models that describe the boundary effects 
in spin chains. 

In this paper we examined the finite-size effects in quasi 
one-dimensional Sr 2 Cu 03 and SrCu 02 - The main objec¬ 
tive of the present study is to carry out an independent 
verification of the field theory result by Sirker et al . n , 
which quantifies the effect of boundaries on the suscepti¬ 
bility of a spin chain. We also investigated the effective 
spin state of Ni 2+ ion doped in a spin \ chain. This ques¬ 
tion is important because the field theory result by Eggert 
and Affleck 15 predicts that an isolated spin 1 impurity in 
a spin | chain may essentially behaves like a scalar im¬ 
purity (S e ff = 0 ) due to the screening cloud formed by 
its antiferromagnetic coupling to its neighbors. 

The plan of the paper is as follows: we first provide 
the experimental details, including a brief summary of 
the results of our crystal growth experiments which will 
be useful in the later part of the manuscript that deals 
with the data analysis and determination of defect con¬ 
centration in the doped crystals. The susceptibility of 
pristine Sr 2 Cu 03 , containing intrinsic chain breaks aris¬ 
ing due to oxygen off-stoichiometry, is analyzed first as it 
presents the simplest case due to its linear chain geometry 
and also allows for a comparison of our results with the 
previous study by Sirker et al. 11 . The procedure is then 
extended to analyze the susceptibility of pristine SrCu0 2 
(zigzag chains) and the Zn and Ni doped Sr 2 Cu 03 crys¬ 
tals. In the zigzag chains, the role of spin frustration 
on defect concentration determination is investigated. In 
the Zn doped Sr 2 CuC >3 crystals, the defect concentration 
obtained is discussed in light of the crystal growth exper¬ 
iment results. The magnetization of the Ni-doped crystal 
is fitted by assuming effective Ni 2+ spin to be zero. The 
influence of Ni and Zn doping on the low-temperature 
phase transition in Sr 2 Cu 03 is presented in the last sec¬ 
tion. 


II. EXPERIMENTAL DETAILS 

The investigations were carried out on high-quality sin¬ 
gle crystals grown using a four-mirror infrared image fur¬ 
nace (Crystal System Corporation, Japan). For the sin¬ 
gle crystal growth of Sr 2 Cu0 3 and its doped variants, the 
starting precursors used and their associated purities (in¬ 
dicated in the parenthesis) were SrC 03 (99.99%), CuO 
(99.995%), ZnO (99.99%), and NiO (99.99%) and in the 
preparation of SrCu0 2 the precursors used were SrCC >3 
(99.9%) and CuO (99.99%). All the successfully grown 
crystals were characterized using optical microscopy un¬ 
der polarized light (Carl Zeiss), scanning electron mi¬ 
croscopy (SEM) (Zeiss Ultra Plus), energy dispersive x- 
ray (EDX) analysis (Oxford), and x-ray powder diffrac- 



FIG. 1 . Unit cell of (a) Sr2Cu03 and (b) SrCu 0 2 ; the crys¬ 
tallographic directions a , b , and c are shown. Representative 
images of the crystals are shown alongside. Lower panels show 
(a) the linear and (b) the zigzag chains. 

tion (Bruker D 8 Advance) obtained by grinding the crys¬ 
tal pieces. The crystals were oriented using a Laue cam¬ 
era equipped with microfocal x-ray source (Photonics Sci¬ 
ence). The dc magnetization measurements were carried 
out along the three principal axes of the grown crys¬ 
tals using a vibrating sample magnetometer (VSM) in 
a physical property measurement system (Quantum De¬ 
sign, USA). Prior to magnetic measurements the crystal 
specimens were annealed under Ar-flow for 72 h at 875°C. 


III. CRYSTAL GROWTH AND STRUCTURAL 
DETAILS 

Sr 2 Cu 03 and SrCu0 2 are known to melt non- 
congruently. Therefore, the crystals of these compounds 
and their doped variants were grown using the traveling- 
solvent-floating-zone (TSFZ) technique associated with 
an image furnace. High purity precursor materials (pu¬ 
rity > 99.99%) were used in the fabrication of these 
crystals, which, combined with the fact that TSFZ is 
a crucible free technique, resulted in high purity single 
crystals, which minimized the effect of extrinsic mag¬ 
netic impurities. Representative images of some of the 
grown crystals with mirror-finished cleaved surfaces are 
shown in Fig. 1. This being a crucible-free technique, 
crystals grown using this process are generally of very 
high-purity, free from contamination from the flux or the 
crucible material. The growth details of SrCu0 2 were re¬ 
cently reported by us in ref. 16. The crystals of Sr 2 CuC >3 
and Sr 2 Cu 0 . 99 M 0 . 01 O 3 (M = Zn, Ni) were grown analo¬ 
gously and their details will be reported elsewhere. Here, 
we briefly present the growth conditions and some crucial 
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TABLE I. Lattice parameters of the grown crystals 


Sample 

a 

b 

c 

Cell volume 

Sr 2 Cu0 3 

12.709 

3.914 

3.499 

174.0 

Sr 2 Cu 0 . 99 Zn 0 . 01 O 3 

12.708 

3.913 

3.499 

174.0 

Sr 2 Cu 0 . 99 Ni 0 . 01 O 3 

12.711 

3.905 

3.498 

173.7 

SrCu0 2 

3.574 

16.332 

3.911 

228.3 


observations pertaining to the segregation of dopants in 
these crystals. 

The crystals were grown under a flowing O 2 atmo¬ 
sphere at a typical growth speed of 1 mm/h. The feed 
and the seed rods were rotated in the opposite direc¬ 
tions at a typical speed of 15-20 rpm to ensure tem¬ 
perature/chemical homogeneity of the floating zone. In 
the growth experiments of the doped crystals we inves¬ 
tigated the composition of the ’’final” floating zone, i.e., 
the floating zone at the end of an experiment that had 
run successfully for several days yielding about 7-8 cm of 
grown crystal. The composition of the final floating zone 
was obtained as follows. The final zone was solidified by 
quickly decreasing the furnace power. The solidified zone 
was sectioned along its length (i.e., parallel to the crystal 
length) and polished to obtain mirror finished surfaces. 
These surfaces were examined using SEM and the aver¬ 
age composition was determined using the EDX probe by 
performing measurements over eight distinct areas span¬ 
ning the entire zone area. In the growth experiment of 
Sr 2 Cu 0 . 99 Zn 0 . 01 O 3 , the final zone composition contained 
11.1 ± 3.5 at. % of Zn, which is significantly larger than 
expected. On the other hand, no trace of Ni accumula¬ 
tion was detected in the corresponding Ni-doped growth 
experiment. The accumulation of solute in the floating 
zone typically arises when the segregation coefficient of 
solute, defined as k = c c lc m where c c / m is the concen¬ 
tration of solute in the crystal/melt, happens to be less 
than 1. In such cases, the concentration of the solute 
in the crystal varies along its length and can be con¬ 
siderably smaller than the nominal value 17 . From our 
observation of Zn accumulation in the floating zone it is 
inferred that the Zn concentration in our grown crystal 
is smaller than the nominal 1 % in the starting material. 
Similar behavior has been reported in the floating-zone 
growth of Zn doped CuGeC >3 crystal 18 . The resolution 
limit of the EDX probe (about 1-2 at. %) does not al¬ 
low us to estimate the doping concentration in the grown 
crystals accurately. However, as shown in the later part 
of our manuscript, chain-break concentration obtained 
from susceptibility analysis also points at a significantly 
smaller Zn doping than nominal 1% in the grown crystal. 

The power x-ray diffraction shows that all the grown 
crystals were phase pure crystallizing with the or¬ 
thorhombic space group Immm for Sr 2 Cu 03 and its 
doped variants, and Cmcm for SrCuCD- The lattice pa¬ 
rameters of all the investigated crystals are summarized 
in Table I. 

The crystal structures are shown in Fig. 1. Sr 2 CuC >3 



FIG. 2. Temperature variation of magnetic susceptibility y: 
(a) Sr 2 Cu 03 and (b) SrCuCL. In both cases data shown are 
recorded by applying field perpendicular to the CuC >2 plaque- 
ttes. Inset: low temperature susceptibilities along crystallo¬ 
graphic a, b, and c axes. 


consists of linear -Cu-O-Cu- chains running parallel to 
the crystallographic b axis. In SrCu 02 , the chain ge¬ 
ometry is zigzag, which can be viewed as a double-chain 
parallel to the c axis that are coupled via a 90° -Cu-O- 
Cu- bond. The intrachain antiferromagnetic interaction 
(J) mediated via the linear -Cu-O-Cu- bond in both these 
compounds is estimated to be of the order of 2000 K. In 
the zigzag chains, the additional intrachain ferromagnetic 
J f via 90° -Cu-O-Cu- bond, was estimated to be about 
0.1J 19 . While the exact role of on magnetic and ther¬ 
modynamic properties is not very clear, probably due 
to its relatively weak strength compared to J, it seems 
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to have no significant influence on the broad features of 
the spin excitation spectrum which shows the expected 
two-spinon gapless continuum at T = 10 K 20 . These 
observations combined with the fact that the interchain 
interaction (J’) is rather weak in both these compounds 
make them excellent realizations of the ID S = | HAF 
model. 


IV. RESULTS AND DISCUSSIONS 

Figure 2 shows the low-temperature magnetic suscep¬ 
tibilities of the pristine Sr 2 Cu 03 and SrCu0 2 crystals. 
The measurements were carried out under an applied 
magnetic field of H = 10 kOe. The data show good agree¬ 
ment with a previous report on well-annealed single crys¬ 
tals due to Motoyama et a/. 14 . The small temperature- 
independent anisotropy above T « 20 K is also in agree¬ 
ment with the previous report and is attributed to the 
Van Vleck paramagnetism 14 . At lower temperatures the 
presence of interchain interaction alters this behavior 
considerably. In Si' 2 Cu 03 , for instance, a distinct change 
of slope of x(T) is observed upon cooling below T ss 5 
K, which coincides with the onset of long-range spin or¬ 
dering in this compound 21 . In SrCu 02 , the absence of 
a similar behavior is in agreement with previous suscep¬ 
tibility and specific heat investigations 22 . Due to inter¬ 
chain interactions, the requirement of good one dimen¬ 
sionality breaks down at low temperature. Judging from 
the susceptibility anisotropy which becomes temperature 
dependent upon cooling below T = 20 K, we expect the 
interchain interactions to be important below this tem¬ 
perature. Since these interactions are not included in our 
theoretical model, therefore, in analysing the susceptibil¬ 
ity for defect concentration analysis we considered only 
the data above T = 20 K. 

As seen in the main panels of Fig. 2, upon cool¬ 
ing below room temperature, susceptibilities of both the 
compounds decrease gradually as expected since at room 
temperature the sample is already cooled far below the 
temperature T = 0.675 J/ks where x(T) is expected 
to show a broad susceptibility peak 2,4 . However, in 


the temperature range below about 100 K or so, the 
susceptibilities of both the compounds show a signifi¬ 
cant ” Curie-tail” arising primarily from the chain breaks 
present due to slight oxygen excess in the crystals 12,14 . 
In the presence of chain breaks one can decompose the 
experimentally measured susceptibility as follows: x = 
Xvv+Xcore+XCurie+Xchain, where Xvv is the Van-Vleck 
paramagnetic contribution and Xcore is the diamagnetic 
contribution of the ionic cores: both these terms are small 
and nearly temperature independent. Henceforth, we will 
use the symbol xo to represent the algebraic sum of these 
two terms; XCurie is the contribution arising due to the 
chain breaks as discussed above and Xchain is the sus¬ 
ceptibility of an uniform endless S = \ HAF chain in 
the absence of defects. Statistically, due to chain-breaks 
there will be as many odd-length segments as even. To 
determine the chain break concentration (p), one can as¬ 
sume that at sufficiently low temperatures, the contri¬ 
bution of the even-length segments dies out and that the 
odd-length segments, having concentration p/ 2 , lock into 
a spin-doublet ground state (Sz = ±|) due the nearest 
neighbor AFM interaction 12 . By fitting the low tem¬ 
perature data using the Curie-Weiss (CW) law one can 
therefore get the value of p. However, this naive approach 
leads to a considerable underestimation of the impurity 
concentration 11 . Moreover , in real systems, a satisfac¬ 
tory fit to the data can only be obtained by including 
the Curie temperature in the fitting expression. How¬ 
ever, the physical significance of the Curie-temperature 
remains ambiguous 12,23 . Fujimoto and Eggert pointed 
out that in the presence of chain breaks, besides the 1/T 
contribution due to odd-length segments, one should also 
take into account the boundary contribution arising from 
the staggered moment at the free ends of a segment which 
also corresponds to a Curie-like behavior with logarith¬ 
mic corrections 10 . Following their work, Sirker et al. ob¬ 
tained a formula for fitting the experimental susceptibil¬ 
ity of a S = \ HAF chain in the presence of chain breaks 
including the boundary corrections 24 . The equation for 
the average susceptibility obtained by them is as follows: 


Xp= (2 -p ) ^ 1 “ - ^ J/T ) + C 1 - P) J/T (( 1 - P + Y^ Xchain + PXB ) (!) 


Here, p is the chain-break concentration and the fac¬ 
tor (1 — p) J / T facilitates the crossover between high- and 
low-tenrperature limits, such that, for a given value of 
p and J, Xp tends to a value ^p^-p) f° r 7 “C P j he., 
at very low temperatures. In this limit, Xp varies as 1/T 
with an effective paramagnetic Curie constant ~ | 

which amounts to an effective spin / paramagnetic con¬ 
centration Ns = This is the same old result obtained 


naively by considering only the odd-length segments con¬ 
tributing to the Curie tail. At the other extreme, i.e., 
■j p (which is easily reached at temperatures well 
below J for small p), the first term diminishes quickly 
to zero; however, the Curie behavior persists due to the 
presence of term pxn = p/[12TZn(2.9J/T)] which cor¬ 
responds to the boundary susceptibility 10 . In this case 
Ns has a logarithmic temperature dependence given by 
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p/[3Tln(2.9J/T)\. At any intermediate temperature the 
effective paramagnetic concentration can be calculated 
using Eq. 7.3 of Ref. 24 . It should be noted that the fac¬ 
tor Na ^ b >1b (here Na is the Avogadro’s constant, g repre¬ 
sents the Lande-g factor, ps is the Bohr magneton, and 
ks is the Boltzmann constant) is not explicitly shown in 
the equation above for Xp ■ I n previous studies, the sus¬ 
ceptibility of these compounds is analyzed by taking a g 
value close to 2 11,12,25 , since the spin associated with the 
Cu 2+ ion behaves like an isotropic Heisenberg spin, which 
is also reflected in the magnetic susceptibilities that show 
a nearly temperature-independent anisotropy. The ESR 
experiments 26,27 also suggest that the value of g for these 
compounds is not significantly different from 2. In view 
of this, in analyzing the susceptibility data using Eq. (1) 
we fixed the g value at 2. Slight departures from this 
value will not change our results significantly. 

Xchain in Eq. (1) is replaced by the fit function ob¬ 
tained by Johnston et al. [Ref. 28 , fit 1, Eq. (50)] 


X.chain — 


4 T 


i + E 
i + E 


q N n 1 
n—1 (T / J) n 
r 

n—1 ( T/J) n _ 


( 2 ) 


The values of N™s and D ra ’s are taken from ref. 28 . The 
expression yields highly accurate values of Xchain in the 
temperature range 0.01J < T < 5J). Since, J for the com¬ 
pounds under consideration is of the order of 2000 K, this 
expression is quite appropriate in the temperature range 
over which the data has been analyzed in the present 
work. It is worth mentioning here that the numerical re¬ 
sult for Xchain due to BF 2 which is frequently used in the 
literature to fit the experimental data of spin chains is 
accurate only at high temperatures T > OAJ/ks which 
makes its applicability inappropriate in the present case. 
This non-applicability of the BF model in fitting the sus¬ 
ceptibilities of these compounds has been demonstrated 
by Motoyama et al. 14 . 

Fitting procedure: In Eq. (1), there are three fit¬ 
ting parameters: (1) the defect concentration p, (2) the 
intrachain exchange J and (3) the temperature indepen¬ 
dent contribution to the susceptibility, Xo {Xp = X~ Xo)- 
Fitting the experimental data to Eq. (1) can give rise 
to arbitrary results if the parameters are not initialized 
judiciously. This problem is further compounded by the 
fact that in the temperature range of our measurements 
the susceptibility values are rather small (« 5 x 10 -5 
emu/mol.Oe). In order to avoid these difficulties, we de¬ 
termined these parameters in steps by selecting suitable 
temperature ranges where Eq. (1) can be further simpli¬ 
fied and the number of parameters reduced. 

In the first step, we considered a temperature range 
over which the measured susceptibility is largely ’p’ in¬ 
dependent which allowed for the estimation of Xo and J. 
For example, in the high temperature range 250 K < T < 
300A" of our measurements the contribution of the first 
term (« p/8T, where p < 0.01 ) is negligible. This leaves 
us with only the Xchain and the Xb terms to deal with 
in Eq. 1. We further note that the relative magnitude 



T (K) 


FIG. 3. Magnetic susceptibility (x — Xo)T of Sr 2 CuC> 3 . Solid 
lines represent the theoretical fit using Eq. (1) (see text for 
details). Data along the b and the c axes are shifted by mul¬ 
tiples of 3 x 10 -3 . 


of pxb is very small compared to (1 - p + pJ/T)xchain- 
For instance, if we assume p to have a value of 0.005 and 
take J to be around 2000 K then near T = 300A” the 
value of the term pxb ~ 7 x 10~ 7 emu/mol.Oe and that 
of (1— p+pJ/T)xchain ~ 9x 10 _5 emu/mol.Oe, where the 
value of Xchain itself is ss 8.7 x 10 -5 emu/mol.Oe which is 
almost 97 % of (1— p+pJ/T)x c hain- Therefore, by taking 
the initial values of J and xo from ref. 14 , we fitted the 
measured susceptibility in the above temperature range 
simply by: x = Xo + Xchain- This gave us a fairly good 
estimate on the values of xo an d J- In the next step, we 
used these values to fit the data from 20 to 300 K using 
Eq. (1) by allowing p to vary. In the final step, all three 
parameters were varied around their values obtained in 
the previous steps to estimate the size of the error bars. 

Susceptibility of Sr 2 CuOz crystal : In Fig. 3 we 
show the temperature variation of quantity (x — Xo)T 
along with the fitted curves. The data and the fitted 
curves are shown on a log-log plot to emphasize the low- 
temperature behavior. A satisfactory agreement of the 
experimental data with the fitted curve is observed over 
the entire fitting range. However, when extrapolated to 
lower temperatures the fitted curve tends to saturate to 
the value ~ §, whereas the experimental data 

continue to decrease due to the intrachain exchange not 
included in our fitting equation. The p values along the 
three crystallographic axes are listed in table II. The av¬ 
erage p value came out to be to be 0.005 per Cu site. 
This value is consistent with a value of 0.006 per Cu site 
reported by Sirker et al. for their unannealed crystal. An 
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TABLE II. Values of intra-chain exchange J, the temperature independent contribution xo and defect concentration p obtained 
by fitting the experimental susceptibility using Eq. (1). For comparison, the low-temperature data were also fitted to the 
Curie-Weiss equation. p_CW and dp are the fitting parameters of this fit (see text for details). 


Parameters 



J 

(k B K) 

Xo 

(10 -5 emu/mol.Oe) 

P 

(per Cu site) 

Pew 

(per Cu site) 

dp 

(K) 


(H 

II a) 


-8.0T0.2 

0.005 

0.001 

2.9 

Sr 2 Cu0 3 

(H 

II b) 

2100 T200 

-7.7±0.2 

0.005 

0.001 

-1.5 


(H 

II c) 


-3.6T0.2 

0.004 

0.001 

-4.8 


(H 

II a) 


1.3±0.2 

0.007 

0.003 

-2.7 

SrCu0 2 

(H 

II b) 

1930±200 

-4.3T0.2 

0.008 

0.003 

-0.7 


(H 

II c) 


-4.3T0.2 

0.008 

0.003 

-3.5 

Sr 2 Cuo.99Zno.oi0 3 

Sr 2 Cuo.99Nio.oi0 3 

(H 

(H 

II c) 

II c) 

2100 T200 

-3.1T0.2 

-3.5±0.2 

0.011 

0.018 

0.004 

0.006 

-1.5 

7.1 


average J value of 2100^200 K and \o values (in the units 
of 10- 5 emu/mol.Oe) of -8.0 ± 0.2(xg), -7.7 ± 0.2(xo) 
and —3.6 ± 0.2(x§), where the superscript on xo denotes 
the crystal axis along which the external magnetic field is 
applied are well within the expected range based on the 
previous studies. For example, Motoyama et al. obtained 
a value of J = 22 OOI 200 by analyzing the susceptibility 
data of annealed crystals up to temperatures as high as 
T = 800A' 14 . Here, it should be mentioned that some¬ 
what higher values of J were reported from the optical 
absorption (J = 2850 K) 29 and the low-temperature spe¬ 
cific heat data (J = 2500 K) 30 , respectively. However, 
based on a Quantum Monte Carlo study it has been sug¬ 
gested that the value of J estimated from y(T) data may 
be reduced due to phonon induced fluctuations of J not 
included in fitting the x(T) data 31 . 

In Table II we also list the chain-break concentra¬ 
tion obtained by fitting the susceptibility data between 
T = 20K and 40 K using the CW law. To fit the CW 
equation, we subtracted the Xo and Xchain contributions 
from the raw data. Xchain is calculated using the J value 
listed in Table II. An average chain-break concentration 
of 0.0012 per Cu site obtained from the CW fit agrees 
well with that reported previously 14,23 , where a similar 
method has been used. However, as expected, this value 
is considerably smaller than that obtained using a more 
rigorous fitting procedure where the boundary contribu¬ 
tions are also taken into account. 

As noted earlier, the effective paramagnetic concentra¬ 
tion (Ns) in our model is a temperature dependent quan¬ 
tity that shows a non-monotonic behavior, decreasing 
rapidly upon increasing the temperature (from a value of 
Ns = vl 2 at T = 0) before increasing again at high tem¬ 
peratures where it goes as p[3ln(2.9J/T) } . At any inter¬ 
mediate temperature T, Ns can be estimated from % p (T) 
by subtracting from it the chain susceptibility Xchain(T) 
and then by equating the temperature times this differ¬ 
ence to an effective temperature dependent Curie con¬ 
stant C(T) (see Eq. 7.3 of Ref. 24 ). Using C(T), the 
value of Ns can be readily calculated (see, for example, 
Ref. 32), using C = (N s N J 4g 2 ^^S(S+l))/3k s . In table 
III, we list the values of Ns for our Sr 2 Cu 03 crystals at 


TABLE III. Effective paramagnetic impurities in the pristine 
Sr 2 Cu0 3 . 


T (K) 

Ns (10 4 per Cu site) 

20 

6.7 

40 

5.1 

80 

4.5 

300 

5.3 


temperatures T = 20,40,80, 200 and 300AL In Fig. 4, 
we compare the calculated M(H) consisting of a Brillouin 
term accounting for the S = b paramagnetic impurities 
(numbering Ns per Cu mol) and the terms linear in H 
consisting of the contribution due to xo and Xchain- The 
calculated curves show a good overall agreement with the 
experimental data at all temperatures and upto the high¬ 
est applied field. The experimental data, however, show 
a weak metamagneticlike transition which is evident from 
a peak in the dM/dH plots shown in the insets. The peak 
position shifts to higher field values with increasing tem¬ 
perature, falling beyond the measurement range of our 
experimental setup at T = 300 K, where M(H) shows a 
linear behavior upto H = 80 kOe, in good agreement with 
the calculated curve. We shall continue our discussion of 
this field dependence below. Here, to conclude this sec¬ 
tion, we note that the values of the fitting parameters 
obtained from our x(T) analysis provide a good overall 
description of the isothermal magnetization at different 
temperatures that span a fairly wide range thereby rein¬ 
forcing confidence in these values. In the future, it will 
be interesting to investigate how this weak metamagnet¬ 
iclike feature seen in M(H) affects the quality of x(T) fit 
at low-temperatures. 

Field dependence of magnetization in S^CuO?,: 

To probe the observed metamagneticlike behavior (Fig. 
4), we measured the temperature variation of magne¬ 
tization under constant applied magnetic fields of in¬ 
creasing strength. The resulting data, plotted as M/H 
against T, is shown in the panel (e) of Fig. 4. Un¬ 
der low fields (H = 10 kOe), the data show a steep 
increase below T = 5K which coincides with the mag- 
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FIG. 4. (a)-(d) Isothermal magnetization (M) of Sr 2 CuC >3 shown as function of applied magnetic field (H). Solid curves are 
calculated M(H) (see text for details). Inset: derivative (dM/dH) as a function of H. Lines are guide to the eye; (e) temperature 
dependence of M/H under various applied fields. In each panel data are shown for field applied perpendicular to the CuC >2 
plaquettes 


netic ordering temperature. At higher fields ( H ss 20 kOe 
and higher), a broad maximum above the ordering tem¬ 
peratures is observed, whose position shifts to higher 
temperatures with increasing field strength. Given that 
the intrachain coupling is around 2000 K, the observed 
field dependence of the chain magnetization under rela¬ 
tively small applied magnetic fields appears surprising. 
This behavior is confirmed by performing similar mea¬ 
surements along the other two crystallographic direc¬ 
tions (not shown here). While it is difficult to pin down 
the origin of this behavior, it is interesting to note that 
63 Cu NMR spectra recorded on a high-quality crystal 
of Sr 2 CuC >3 also exhibit anomalous features in the same 
field-temperature range 33 . The field sweep NMR spec¬ 
tra develop a broad background with sharp edges at low 
temperatures, which was theoretically predicted as aris¬ 
ing due to external magnetic field induced local staggered 
magnetization near the chain breaks 34 . However, two ad¬ 
ditional features in the NMR spectra, not expected from 
the theory of Ref. 34 , were also observed, viz., splitting 
of the main NMR peak and appearance of shoulders on 
both sides of the main peak. In one of the theoretical 
studies, these features were attributed to ’’mobile bond 
defects,” which can arise due to spin-lattice coupling and 
result in a ’’local alternating magnetization” 35 . In Ref. 
36, however, the role of interchain exchange in giving 
rise to these additional NMR features has been scruti¬ 
nized. In the future, it will be interesting to investigate 
if there is any correspondence between the anomalous 


NMR line-shape and magnetization behavior observed in 
our experiments. 

Susceptibility of SrCuO '2 crystals: Next, we fitted 
the susceptibility of SrCuC> 2 . The fitted curve is shown 
in Fig. 5. A satisfactory fit to the data is obtained down 
to nearly T = 30 K. An average J value of 1930l2oo^ 
and xo = 1.3±0.2( X g), -4.3±0.2(*&), -4.3±0.2(x§) (in 
the units of 10 5 emu/mol.Oe) obtained from the fitting 
analysis are in good agreement with Motoyama et al. 14 . 
The average defect concentration turned out to be 0.008 
per Cu site which is higher than that in Sr 2 Cu 03 which 
was grown and annealed under similar conditions but 
prepared using starting precursors with slightly higher 
purity. 

In Fig. 6 we compare the susceptibilities of the two 
compounds by plotting (y — y 0 )J against T/J where 
the values of yo and J for the two compounds are used 
from table II. In this representation, the susceptibility 
is J-independent. On the same plot we also show the 
variation of XchainJ 6 - Near room temperature the two 
data sets collapse on the calculated Xchain as one would 
have expected. However, upon lowering the temperature 
below T ss 0.1J the measured susceptibilities start to 
deviate from the ideal chain behavior: increasing with 
decreasing temperature due to chain breaks. This in¬ 
crease is more pronounced in SrCu 02 than in Si^CuOs, 
in agreement with the higher p value for the former. The 
double chains of SrCuC >2 are characterized by the pres¬ 
ence of an additional ferromagnetic frustrated exchange 
















FIG. 5. Magnetic susceptibility (x — xo )T of SrCu02. Solid 
lines represent the theoretical fit using Eq. 1 (see text for 
details). Data along b and c axes are shifted by multiples of 
3 x 10“ 3 

Jp which is estimated to about 0.1J 19 . While the dom¬ 
inant antiferromagnetic exchange (J) tends to lower the 
chain susceptibility upon lowering the temperature, Jp 
has just the opposite effect as it tends to disrupt any 
tendency towards antiferromagnetic ordering. Therefore, 
the chain susceptibility Xchain dr SrCuC >2 is expected to 
show a slower rate of decrease upon cooling than for a 
linear chain as was used in Eq. (1) for the estimation 
of p. Also, the presence of coupling in the double 
chains itself implies that a defect cannot be considered 
as a chain break at low temperatures. Therefore, at tem¬ 
peratures below the rung coupling Jj?, Eq. (1) has only a 
limited validity. The p value obtained here for undoped 
SrCuC >2 is, therefore, possibly slightly overestimated. A 
more direct manifestation of evidently is to suppress 
the long-range spin ordering which has been studied pre¬ 
viously using y SR and neutron studies 22 . It is worth 
mentioning here that though the next-nearest neighbor 
interaction J nnn along the chain length is estimated to 
be around 140 K 3 ', Eq. (1) is shown to be applicable 
as long as T> p 2 J n nn 24 - Thus J nnn is not expected to 
have any significant influence either on the susceptibility 
or on the determination of p. 

Susceptibility of Sr2 Cuo.ggZnoxn O3: We investi¬ 
gated the behavior of S 12 CUO 3 crystal doped with non¬ 
magnetic Zn (nominal doping level ~ 1%). As mentioned 
in the crystal growth section, by analyzing the composi¬ 
tion of the frozen-in floating zone we had inferred that 
the Zn concentration in the grown crystal should be sig¬ 
nificantly smaller than the nominal value. To test this 
assertion, we analyzed the susceptibility data of the Zn- 



FIG. 6. Normalized susceptibilities of Sr 2 CuOs and SrCu02- 
In both cases data shown are for field applied along the 
longest crystallographic axis. The calculated chain suscep¬ 
tibility Xchain J is also shown (solid line). The dashed lines 
represent calculated \ V J for p = 0.005 and 0.008 per Cu site 
(see text for details). 


doped crystals using the same procedure as used for the 
pristine crystals. The results are shown in Fig. 7. A 
reasonably good fit to the experimental data is observed 
down to nearly T = 40 K. The value of the exchange inte¬ 
gral did not show any appreciable decrease upon doping, 
which is not surprising given the dilute levels of dop¬ 
ing. Similarly, %o also did not show any significant vari¬ 
ation with respect to the pristine compound (Table II). 
The chain-break concentration was found to be around 
0.011(2) per Cu site. Taking a cue from the work of 
Sirker et ai, where the chain-break concentration in the 
doped crystals exceeded the nominal Pd concentration by 
an amount roughly equal to the concentration of intrin¬ 
sic chain breaks in their pristine compound, we conclude 
that the concentration of chain breaks induced by Zn 
doping in our crystal is around ~ 0.006, the remaining ~ 
0.005 being due to the intrinsic defects as inferred from 
the analysis of the pristine compound. This smaller Zn 
concentration of 0 . 6 % obtained from the susceptibility 
analysis is in agreement with the crystal growth results. 

Susceptibility of Sr 2 Cu 0 . 99 IVi 0.01 O 3 : Recently, it 
was reported that upon 1 % Ni doping in SrCuC >2 a dra¬ 
matic change of the low-energy spin excitations takes 
place. Briefly, what had been a gapless two-spinon 
excitation spectrum for the pristine compound became 
gapped upon Ni doping 38 . The value of the spin pseudo¬ 
gap was found to be around ~ 8 meV. In our own pre- 
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FIG. 7 . Magnetic susceptibility (x ~ Xo)T of 

Sr2Cu0.99Zn0.01O3. Solid curve in the main panel rep¬ 
resents the best theoretical fit using Eq. ( 1 ) 


liminary experiments using inelastic neutron scattering 
we found a similar gap in Ni-doped Sr 2 Cu 03 . However, 
no such gap was found in the Zn 2+ doped crystals 39 . Eg- 
gert and Affleck 15 investigated the behavior of an isolated 
magnetic impurity substituted in a half-odd-integer-spin 
HAF chain theoretically. In particular, they considered 
an isolated S = 1 impurity (e.g., Ni 2+ ) in the chain and 
showed that in the limit of infinite chain length Ni 2+ 
moment either remains decoupled from the chain or is 
completely screened by the antiferromagnetically coupled 
neighboring spins. In the latter case, the screened Ni 2+ 
moment essentially breaks the chain at the impurity site. 

Therefore, if the spin of Ni 2+ ion is indeed screened we 
should be able to get the Ni concentration in our crys¬ 
tal using Eq. (1). The measured susceptibility of the 
Ni-doped crystal was analyzed using Eq. (1). In Fig. 
8 we show the calculated susceptibility for p = 0.018(1) 
obtained by fitting the experimental data using Eq. (1) 
as done for the pristine crystal. Within the accuracy of 
our method, the fitted values of J and xo are the same as 
for the pristine crystal (see, Table II). The defect concen¬ 
tration estimated here is higher than the nominal value 
of 0.01. However, if one takes into account the intrin¬ 
sic defect concentration due to oxygen off-stoichiometry, 
deduced from the analysis of the pristine compound pre¬ 
pared and annealed under the identical conditions, we get 
a value of around 0.013(1), which is indeed close to the 
nominal concentration. We also compared the measured 
M(H) at T = 20 K with that calculated for two differ¬ 
ent spin states: S = 1 and S = 0 for Ni concentration 
of 1 %. The comparison is shown in the inset of Fig. 8. 


FIG. 8. Magnetic susceptibility (x — Xo)T of 

Sr2Cu0.99Ni0.01O3. Solid curve in the main panel rep¬ 
resents the best theoretical fit using Eq. ( 1 ). Inset: 
isothermal magnetization at T = 20 K. Solid curves are 
calculated M(H) for S = 0 and S = 1 cases (see text for 
details). 

The calculated magnetization in the S = 1 case is clearly 
very large and deviates considerably from the experimen¬ 
tal data. On the other hand, in the S = 0 case a close 
match with the experimental data is observed. There is, 
therefore, an overall agreement between the experimental 
data and the calculated curve that proves the point that 
Ni 2 + spins are indeed screened. This is also in agree¬ 
ment with the analysis reported in ref. 25 on the powder 
samples of Ni doped Sr 2 Cu0 3 , where it was shown that 
the Curie tail at low temperatures is considerably smaller 
than expected from unscreened S = 1 Ni moments. 

It should be pointed out that the fit-quality for the 
Ni-doped case using Eq. 1 is not as good as for the pris¬ 
tine crystal, particularly below T = 100 K. The reason 
why the fitting Eq. 1 did not work so well is proba¬ 
bly related to the fact that the ground state of Ni-doped 
crystal is gapped which makes the Xchain used in Eq.(l) 
not suitable at temperatures below which the gap opens. 
In analogy with SrCu0 2 , the gap value is expected to 
be around 8 meV (~ 80 K) for a similar Ni concentra¬ 
tion, which roughly coincides with the temperature below 
which the fit quality deteriorates. To conclude, our anal¬ 
ysis of the magnetization data supports the screening of 
the Ni 2+ moments. 

It is intriguing as to why Ni doping opens up a gap 
while the Zn doping apparently does not. An important 
difference between these two impurity types is that unlike 
the chain breaks due to Zn 2_t_ impurity, where the spins 
on the either side of the impurity may still couple weakly 
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FIG. 9. Cartoon showing screening of the Ni spin (S = 1) 
coupled antiferromagnetically to it neighbors in a spin-| an¬ 
tiferromagnetic chain. 

through the next-nearest-neighbor interaction, the Ni 2+ 
impurity results in a spin compensated cloud around the 
impurity site, which isolates the chain segments (Fig. 9). 
In such a scenario, whose validity needs to be established 
through further work, the spinons will be confined over 
finite-length sections of the chain which should open up 
a gap in the low-energy spectrum with the gap value 
scaling inversely with the average length of the segments 
(L ~ 1). Another possible reason for not having found a 
gap in the Zn 2+ crystals might be related to the fact that 
Zn 2+ does not readily replace Cu 2+ in a crystal grown 
from the melt; hence the quantity of Zn doped in the 
chains may not be enough to open up an observable gap. 
Recently, spin gap has also been reported in the crys¬ 
tals of Sr 1 .9Cao.iCu0 3 40 and Sr 0 .gCao.iCu0 2 41 . Since 
Ca 2+ replaces Sr 2+ , the -Cu-O-Cu- chains in these crys¬ 
tals should remain unsegmented. It has been argued that 
the ionic-size mismatch of Ca 2+ and Sr 2+ gives rise to a 
local bond disorder of the intrachain exchange coupling J 
which is probably responsible for the observed spin gap. 
The appearance of spingap is, therefore, not always as¬ 
sociated with the segmentation of the chains and sub¬ 
stantial bond disorder can also cause the ground state to 
be gapped. In the light of these recent results, it is fair 
to admit that the chain break picture is probably too 
naive. Even the notion of screening, strictly speaking, 
is valid at T = 0. At non-zero temperatures, the sus¬ 
ceptibility will be affected by the screening cloud around 
the impurity spin, which extends beyond nearest neigh¬ 
bors making the applicability of Eq. (1) in such scenarios 
rather restrictive. 

Magnetic ordering in doped SrzCuO?,: Here we 
briefly discuss the influence of doping on the magnetic 
ordering temperature. Previously, the weak magnetic or¬ 
dering in Sr 2 Cu0 3 at T = 5.4 K was reported based on 
neutron and p,SR studies 23 . From neutron studies the 
ordered Cu moment is reported to be < 0.06/is 23 . A 
small but clearly discernible peak near T ss 5.5 K in 
our susceptibility data measured under an applied field 
of 1 kOe (Fig. 10) indicates a transition to an ordered 
state. The antiferromagnetic ordering of the Cu 2+ mo¬ 
ments is also reflected in the isothermal magnetization 
of the pristine sample at T = 2 K for H || b and H || 
c (i.e., parallel to the chains and perpendicular to Cu0 2 



FIG. 10. (a) Low temperature susceptibilities \ of the pris¬ 
tine, Zn and Ni doped Sr 2 CuC >3 crystals. Measurements were 
carried out under an applied field H = 1 kOe along the chain 
direction. Inset: isothermal magnetization M(H) of the pris¬ 
tine crystal along crystrallographic b and c axes. The arrow 
indicates the field induced transition in the H || b measure¬ 
ments. The solid line is the calculated M(H) at T = 2 K (see 
text for details), (b) Specific heat C p /T is shown for the three 
compounds. The lines are guides to the eye. Arrows show the 
magnetic transitions. 


plaquettes, respectively). The data are shown as an inset 
in the panel (a) of Fig. 10. In the ordered ground state, 
M(H) exhibits a sharp field induced transition near H = 
5 kOe along the chain direction (H || b), while along H 
|| c it shows only a monotonic increase. The calculated 
M(H) for N s 1.8 xlO -3 per Cu site, obtained using the 
p value listed in table. II, is compared with the experi- 
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mental data. With increasing field values, the calculated 
curve increases far more rapidly than the experimental 
data. This was expected, since at T = 2 K the spins are 
in an antiferromagnetically ordered state. Upon doping 
Zn the transition temperature is suppressed to T = 4 
K. To confirm the bulk nature of these transitions, spe¬ 
cific heat measurements were performed where anomalies 
were seen at the same temperatures as the susceptibili¬ 
ties (shown in the inset of Fig. 10) where molar specific 
heat divided by temperature ( Cp/T ) is shown as a func¬ 
tion of T. The detection of weak magnetic ordering in 
our bulk measurements is an indication of the high qual¬ 
ity of the single crystals used in the present work. The 
transition temperature upon Pd doping in Si' 2 Cu 03 is 
shown to decrease 23 in agreement with a quantum me¬ 
chanical calculation 42 . The observed suppression of the 
transition temperature upon Zn doping is in a good qual¬ 
itative agreement with this theory. The effect of Ni dop¬ 
ing on the magnetic ordering temperature, however, is 
quite dramatic as no signs of transition to an ordered 
state either in the susceptibility or the specific heat are 
observed down to the lowest measurement temperature 
of T = 2 K in our experiments. According to the theory 
of ref. 42 , the transition temperature under 1 % of scalar 
impurity, should have been suppressed to about 0.5Tjv 
(see, Fig. 7 of ref. 23 ). However, the absence of any sign of 
ordering in our Ni-doped crystal probably indicates that 
opening of a sizable spin gap has ousted the magnetic 
ordering completely or has suppressed it to temperatures 
far below T = 2 K. It should be noticed that the specific 
heat of Ni doped crystal at low temperatures is slightly 
(about 2 — 3 mJ/molK 2 ) smaller in magnitude compared 
to that of the pristine or Zn doped crystal. This decrease 
is probably a manifestation of the spin pseudogap. A 
complete analysis of the specific heat in the absence of 
a lattice template and without the knowledge of contri¬ 
bution arising from the free spins at the chain ends is 
difficult at this stage. 


V. CONCLUSION 

In conclusion, we studied the finite-size effects in the 
quasi-one-dimensional spin S = \ Heisenberg antiferro- 
magnets Sr 2 Cu 03 , Sr 2 Cu 0 . 99 M 0 . 01 O 3 (M = Zn and Ni) 
and SrCu0 2 - The investigations were carried out on 
high-quality single-crystalline samples. We experimen¬ 
tally verified the field theory results published by Sirker 
et al. 11 ’ 24 describing the finite-size effects on the mag¬ 
netic susceptibility of a spin | Heisenberg Antiferromag¬ 
netic chain.We also showed that dilute concentration of 
Ni 2+ ions in the chain behave like scalar defects as pre¬ 
dicted by the theory of ref. 15 , which is probably related 
to gapping of the spinon excitation spectrum upon Ni 
doping as shown by Simutis et al. 38 . We also show the 
inadequacy of the model in describing the susceptibility 
in two cases: (i) when the ground state is gapped and (ii) 
in the zigzag chains characterized by the presence of weak 
rung coupling . 38 Our work should motivate further field 
theoretic studies to extend the model in order to remove 
these inadequacies. In future, experimental verification 
of this theory on other spin chain compounds should be 
done. The influence of doping on the low-temperature 
long-range spin ordering temperature was studied. Zn- 
doping suppresses the magnetic ordering in accordance 
with the theory. On the other hand, in the Ni doped 
crystal no signs of magnetic ordering were found down 
to the lowest measurement temperature of T = 2 K. The 
complete ousting of the magnetic ordering or its suppres¬ 
sion far below T = 2 K is probably due to the gapped 
ground state. In the compound SrCu0 2 , consisting of 
zigzag S = ^ chains, the influence of spin frustration 
on the magnetic ordering and the defect concentration 
determined from the susceptibility data is discussed. 
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